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ABSTRACT 

This paper presents transition of the failure mode of a cohesive, spherical body due 
to YORP spin-up. On the assumption that the distribution of materials in the body 
is homogeneous, failed regions first appearing in the body at different spin rates are 
predicted by comparing the yield condition of an elastic stress in the body. It is found 
that as the spin rate increases, the locations of the failed regions move from the 
equatorial surface to the central region. To avoid such failure modes, the body should 
have higher cohesive strength. The results by this model are consistent with those by 
a plastic finite element model. Then, this model and a two-layered-cohesive model first 
proposed by Hirabayashi et al. are used to classify possible evolution and disruption 
of a spherical body. There are three possible pathways to disruption. First, because of 
a strong structure, failure of the central region is dominant and eventually leads to a 
breakup into multiple components. Second, a weak surface and a weak interior make 
the body oblate. Third, a strong internal core prevents the body from failing and only 
allows surface shedding. This implies that observed failure modes may highly depend 
on the internal structure of an asteroid, which could provide crucial information for 
giving constraints on the physical properties. 

Key words: minor planets, asteroids: general - protoplanet ary discs - methods: 
analytical - methods: numerical 


1 INTRODUCTION 

Recent observations have shown that rotational disruption 
may be common in our system. P/2013 R3 has broken into 
multiple components (Jewitt et al. 2014), 62412 has had 
dust tails from its nucleus (Sheppard & Trujillo 2015), and 
P/2012 F5 (GIBBS) has ejected dust and fragments (Drahus 
et al. 2015). They are considered to be candidates that have 
failed due to fast rotation, implying that rotational disrup¬ 
tion may be diverse. 

The YORP effect, which highly depends on the obliq¬ 
uity and shape (Scheeres 2007; Statler 2009; Cotto-Figueroa 
et al. 2015), may be one of the primary factors that accel¬ 
erate/decelerate the spin of an asteroid with a few hundred 
meters in diameter (Rubincam 2000). We have to mention 
that small impacts are also capable of changing the spin 
rate of an asteroid quasi-statically (Marzari et al. 2011). 
Spin up/down by this effect may be more critical to aster- 
oidal evolution in the solar system (Pravec et al. 2006) and 
in the debris disk around a white dwarf (Veras et al. 2014) 
than ever thought. A remarkable aspect of fast rotating as¬ 
teroids is that observed asteroids spinning at spin periods 
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shorter than 3.5 hours are all spheroidal (personal communi¬ 
cation with Patrick Taylor, 2015). For example, 1999 KW4 
Alpha, which is rotating with a spin period of 2.76 hours, 
is a quasi-spherical body with ridges at its equator (Ostro 
et al. 2006). The similar features were also seen for the cases 
of 1950 DA (Busch et al. 2007) and 1994 CC Alpha (Bro- 
zovic et al. 2011). This may tell us that fast rotating, quasi- 
spherical asteroids are fairly common in the solar system. It 
is also important to mention that rotational disruption by 
the YORP effect would lead to complex dynamical evolution 
(Jacobson & Scheeres 2011). 

Such analytical and observational evidence of quasi¬ 
static spin-up by the YORP effect has been giving rise to the 
question about evolution and disruption of quasi-spherical 
objects. Analytical works have contributed to better under¬ 
standing of failure mechanisms of a uniformly rotating ellip¬ 
soid (Dobrovolskis 1982; Davidsson 2001; Holsapple 2001, 
2004; Sharma 2013). Numerical investigations have been ca¬ 
pable of dynamically modeling shape reconfigurations by the 
YORP spin-up (Walsh et al. 2008, 2012; Sanchez & Scheeres 
2012). These studies gave scientific insights into reshaping 
processes of ellipsoidal and quasi-spherical objects affected 
by material friction. Significant progress was made by Walsh 
et al. (2008) who first succeeded in explaining the formation 
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of equatorial ridges on a quasi-spherical body by using a 
hard-sphere discrete element code. Also, Sanchez & Scheeres 
(2012) used a soft-sphere discrete element code to analyze 
the failure modes of cohesionless ellipsoids having different 
friction. An ideal shape of an initially spherical body affected 
by self-gravity and rotation was numerically and analyti¬ 
cally investigated by considering the angle of repose (Harris 
et al. 2009; Scheeres 2015). Regardless of these efforts, the 
reshaping process of a quasi-spherical body is still poorly un¬ 
derstood. Specifically, how does the internal structure con¬ 
tribute to evolution and disruption due to quasi-static spin- 
up? How does the stress field have an effect on failure at 
severe rotation conditions? How does cohesive strength help 
an asteroid remain without structural damage? 

The motivation of this study comes from earlier studies 
about the failure conditions of a quasi-spherical body spin¬ 
ning fast. Assuming a homogeneous structure, Hirabayashi 
& Scheeres (2015) found that 1950 DA may mainly experi¬ 
ence failure of its central region, but should not have surface 
failure at the current spin period, 2.1 hours. This result in¬ 
dicates that landslides may not make equatorial ridges on 
a quasi-spherical asteroid, but inelastic deformation due to 
strong compression along the spin axis is a likely driver for 
it. Hirabayashi et al. (2015) confirmed that at high spin 
rates the central region certainly reaches the yield first and 
alternatively proposed that to have surface failure and mass 
shedding driven by it a spherical body must have a strong 
core. These studies imply that the internal structure plays 
a prime role in evolution and disruption of a quasi-spherical 
object. The issue in their works, however, is that they only 
considered faster spin cases. 

The main objective in this paper is to infer the de¬ 
tailed failure mode of a quasi-spherical object at different 
spin states and to visualize its evolution and disruption. To 
do so we consider an asteroid to be perfectly sphere and an¬ 
alyze its failure modes and conditions by applying a contin¬ 
uum mechanics approach. Theoretical and numerical studies 
are employed here. For the theoretical investigation, ana¬ 
lytical solutions of the stress field by linear elasticity are 
used to predict when and where failure occurs in a spher¬ 
ical body. For the numerical study, a plastic finite element 
model (FEM) developed by Hirabayashi & Scheeres (2015) 
is employed to see the detailed conditions for failure of the 
body. 

It should be noted that the analytical study here di¬ 
rectly applies the technique by Dobrovolskis (1982). He in¬ 
troduced a quantity of shear and pressure called “duress” 
to investigate the internal structures of the Martian moons, 
Phobos and Deimos, and several asteroids. According to his 
study, Phobos should have the maximum shear stress at the 
central region at the current orbital and rotational conhg- 
urations, and could only have thrust faults on its surface. 
He also studied the fate of Phobos, and found that the tidal 
force will have a primary effect on its failure once it reaches 
a distance of 6700 km from the center of Mars. We extend 
his pioneering work to investigate the case in which a spheri¬ 
cal object changes its stress field due to quasi-static spin-up. 
Specifically, we focus on how possible failure modes can shift 
in the body by considering the Drucker-Prager yield crite¬ 
rion, a pressure-shear dependent yield condition. 

This paper is organized as follows. First, using linear 
elasticity, we establish a stress field in a spherical body as 


a function of a spin rate. Second, the yield conditions of 
the surface region and the central region are analytically 
investigated to predict the locations of the failed regions that 
appear first. Third, we use a plastic FEM to evaluate the 
analytical results. Finally, possible evolution and disruption 
scenarios of a spherical body are proposed by taking account 
of a two-layered-cohesive model by Hirabayashi et al. (2015). 


2 MODELING OF A SPHERICAL ASTEROID 

We begin our analysis with modeling of a spherical asteroid 
that experiences the YORP spin-up. A spherical object is 
supposed to have a homogeneous structure. Although the 
model used is the same as Hirabayashi et al. (2015), we do 
not explicitly consider heterogeneity in the body as our focus 
is on determining the critical values of material properties 
with which local regions fail structurally. The density, the 
radius and the gravitational constant are defined as p, R and 
G, respectively. The following discussion will use dimension¬ 
less quantities. Lengths, body forces, spin rates and stress 
tensors are normalized by R, ivpGR, y/TvpG and ixp^GR?, 
respectively. The normalized spin rate and radial distance 
from the center are denoted as uj and Rh, respectively. 

The spherical body is assumed to be rotating along a 
constant spin axis. Here, this spin axis is lined up in the « 
axis. The x axis is then given in an arbitrary direction on the 
equatorial plane, and the y axis is defined to satisfy the or¬ 
thogonal relation with these axes. This Cartesian coordinate 
frame is later denoted as {x,y,z). We also use a spherical 
coordinate system {Rb,8,tp), where 0 is the latitude and ip 
is the longitude. 

Material behavior is assumed to be elastic-perfectly 
plastic - that is, linearly elastic below the yield condition 
and no-material hardening on it. To describe such behavior, 
we use three material properties: Poisson’s ratio u, a fric¬ 
tion angle (p a-nd cohesive strength Y. Although we evaluate 
a stress field based on linear elasticity, since our focus is only 
on the stress. Young’s modulus does not appear in our for¬ 
mulation. Thus, we do not specify this value here. We also 
consider the latitude 9 a free parameter in the theoretical 
investigation. 

In this elastic-perfectly plastic model, only elastic strain 
contributes to stress (Simo & Hughes 1998). Also, at the 
yield condition stress does not change, and thus only plastic 
strain occurs (Chen & Han 1988). On the assumption that 
plastic deformation is small, this fact guarantees that the 
yield condition of an elastic stress predicts when and where 
plastic deformation appears and propagates. To calculate 
an elastic solution for a spherical body, again, we use the 
analytical form derived by Dobrovolskis (1982). 

The present study needs to take into account the fact 
that plastic deformation is a function of loading paths, which 
are time-dependent. However, it takes about million years 
for the YORP effect substantially to change the spin rate 
of an asteroid (Rubincam 2000), implying that the changes 
due to time are negligible. A recent breakthrough about the 
stochastic YORP effect may also provide an interpretation 
that an asteroid would rather change its spin state stochas¬ 
tically than spin up constantly (Statler 2009; Bottke et al. 
2015). This intimates that the internal structure may go 
back and forth between elastic and plastic states. 
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Holsapple (2010) who analyzed plastic deformation of 
a body due to the YORP effect simplified his analysis with 
the following crucial statement, “Clearly the possible anal¬ 
yses are vast in parameter space, and are much too compli¬ 
cated to expect closed-form solutions. A much simpler and 
reasonable approach is to ignore the details of the applica¬ 
tion and timing of those surface forces, and study the case 
where the angular momentum is slowly increased in some 
unspecified but quasi-static way.” We follow his assumption 
to consider the present problem a quasi-static problem. We 
leave analyses about dynamical evolution of the stress field 
in an asteroid as our future work. 

We define the terminology about cohesive strength used 
in this study. The critical cohesive strength is a computed 
threshold that can prevent an asteroid from experiencing 
failure. The “actual” cohesive strength, on the other hand, 
is an assumed quantity that an asteroid may actually have, 
which will be constrained by the critical cohesive strength. 
Also, the spin rate, uj = ~ 1.15, will be used to de¬ 

scribe the condition at which a small particle resting on the 
surface of a sphere has zero net acceleration (Hirabayashi 
et al. 2015). This simply implies that beyond this spin rate, 
materials there experience tension and are shed. 


3 ANALYTICAL INVESTIGATION 
3.1 Yield condition 

Although Dobrovolskis (1982) defined duress to express the 
dependence of materials on pressure and shear, here, we use 
the Drucker-Prager yield criterion, a smooth function that 
takes account of shear and pressure: 

f = ah + - s ^ 0. (1) 


h and J 2 are the stress invariants: 

h = (Ti+a2 + a3, (2) 

J2 = ^{(CTI - CT2)^ -I- (0-2 - 0-3)^ -I- (0-3 - CTl)^}, ( 3 ) 

where Ui (i = 1, 2,3) is the principal stress component, a and 
s, which are functions of (f) and Y, characterize the shape of 
this yield surface. In general, the radius on the deviatoric 
stress plane is determined by choosing appropriate a and s. 
To determine these parameters, Chen & Han (1988) com¬ 
pared the radius of the Drucker-Prager criterion with that 
of the Mohr-Coulomb yield criterion. Here, we use the fact 
that the stress field in a spherical object is always located at 
the compression meridian, i.e., ai — a 2 > as (Hirabayashi 
& Scheeres 2015). Then, a and s are (Chen & Han 1988): 

2 sin (A OYcoscA 

a = —= -, s = —=-. (4) 

%/3(3 — sin(()) \/3(3 — sin (()) 


3.2 Elastic stress solution in a spherical body 

Hirabayashi et al. (2015) used a technique by Dobrovolskis 
(1982) to provide an elastic stress in a spherical body rotat¬ 
ing at spin rate iu. Here, we directly apply their formulation 
to the present problem. Since the body forces acting on the 
body are symmetric along the 2 axis, we consider the stress 
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field at -!/) = 0: 


(IFexx 

= k\{l — Rb) — kuRl sin^ 9, 

( 5 ) 

CTeyy 

= ki {1 — rI) — kisRl cos^ 9— 

kuRb sin^ 9, (6) 

^ezz 

= kg{l — rI) — ki4,Rlcos^ 9, 

( 7 ) 

^ exz 

where 

= ki^Rl cos 9 sin 9. 

( 8 ) 


fci = 

ui'^ 12+i^-5i^^ 

2 3-iv 

5 (l-i/)(7+5v) 

15 1-iv ’ 

kg = 

^'2 3 6^5,y2 

2 3-1^ 

5 (l-i/)(7+5i/) 

15 l-i^ ’ 

II 

CO 

2lj^ 

1 4 1-21/ 

5 (1-i^)(7+52/) 

3-6u-5v^ 

' 15 1-u 

1 4 l-2u 

15 1-1/ 

II 

5 (l-i^)(7+5i^) 


( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 


The indices of ks are chosen to be consistent with those by 
Hirabayashi et al. (2015). 


3.3 Characterization of failure conditions 


The yield condition of the stress field is investigated to char¬ 
acterize when and where failure occurs in the body. Specifi¬ 
cally, we consider the critical cohesive strength, denoted as 
Y*. From the equal condition of Form (1), we obtain Y* as 
a function of lo, (f>, 9, Ri, and v. 

Y* {lj, 0 , 9, Rt, u) = ~ '^A (q7i + ^). ( 13 ) 

6 cos (p 

As discussed in the previous section, since we assume that a 
material is elastic-perfectly plastic, if Y* at a given location 
is higher than the actual cohesive strength, that location 
reaches the yield. Thus, giving the contour curves of Y* 
predicts the failed regions in the body. A negative value 
of Y* indicates that the locations having it do not need 
cohesive strength to avoid the yield. 

We analytically explore how Y* evolves as to changes. 
To do so we compute Y*s at the body’s origin and on the 
surface, later denoted as Y* and Y/, respectively. At the 
origin, where Rb = 0, the stress invariants are: 


h — 2/ci -f kg , 

J 2 = ^(fci-fcg)". 
Then, we obtain Y)*: 


Y*{LU,(j>,u) 


3 — sin 4> Z + 2v 2 
6 cos (j) 7 + 


(14) 

(15) 


^ {3-iy)(-2 + ujh 

-I-Q 

5 1-u 


(16) 


Needless to say, Y* is independent of 9. When Y* = 0, the 
spin rate should be: 


LJ 


*2 

c 


2 

1 I 5 l — 3+2i/ 

\/3q: 3 — u 7+5i/ 


(17) 


This condition is a key parameter that tells us when the 
central region needs cohesive strength to avoid failure. This 
parameter indicates a lower threshold showing when a spher¬ 
ical body becomes oblate. 

On the surface, where Rt ~ 1, on the other hand, the 
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stress invariants are: 

Ii = —(fci3 + fci4) cos^ 0 — 2fci4 sin^ (18) 

J 2 = ^{(fci3 — fci3fci4 + fci4) cos"^ 0 + A:i 4 sin'^ 6» 

+fci4(A:i3 + fci3) sin^ 0cos^ 0}. (19) 

Then, T/ and the spin rate at which T/ = 0, denoted as cOg, 
are obtained as functions of 6 by substituting these equa¬ 
tions into Equation (13); however, since the expression is too 
long to write down here, we only mention that the derivation 
is trivial. 


3.4 Results 

Hereafter, following the fact that the failure conditions are 
not strong functions of Poisson’s ratio (Holsapple 2008; 
Hirabayashi & Scheeres 2014), we fix Poisson’s ratio at 0.25, 
which allows for compressibility of the volume. For <j), we 
consider a range of 25° — 45° for a typical geological mate¬ 
rial (Lambe & Whitman 1969). 

Figure 1 shows the contour plots of W* (the dashed 
curves) and Y* (the solid curves) in a, 9 — co space. Figures 
1(a), 1(b) and 1(c) display the cases of (j) = 25°, 35° and 45°, 
respectively. The red solid and dashed curves provide tuj 
and Wc, respectively. The blue curves show the condition at 
which Y* = Y*, which means that if the spin rate is higher 
(lower), the central region is more (less) sensitive to failure 
than the surface. The fact that the blue line and are close 
to one another shows how rapidly Y* increases. Thus, can 
be considered to a parameter describing when the central 
region starts to fail structurally. For the cases of <j) = 25° 
and 35°, Y* is always positive. The analytical expression of 
y/ shows that it is always positive if 1 — 12a^ > 0, i.e., 
(j) < sin-i(3/5) ~ 36.87°. 

We plot the distributions of Y* on the x — z plane for 
different spin rates and friction angles in Figure 2, which are 
analogous to Figure 4 in Dobrovolskis (1982) except that he 
displayed the maximum shear. The left and right columns 
indicate spin rates of 0.9 and 1.15, respectively. The 0.9 spin 
rate is chosen based on uj*. Figures 2(a) and 2(b), 2(c) and 
2(d), and 2(e) and 2(f) display the cases of friction angles of 
25°, 35°, and 45°, respectively. 

For the case of a friction angle of 25°, = 0.85 and 

Y* > 0 everywhere. At a spin rate of 0.9, if the actual cohe¬ 
sive strength, later simply denoted as Y, is less than 0.02, a 
region enclosing the surface and the central region becomes 
the most sensitive to failure (Figure 2(a)). We note that if 
the spin rate is less than w*, the central region should be 
stable. At a spin rate of 1.15, however, since Y* becomes 
maximum at the central region, if Y is less than, for exam¬ 
ple, 0.175, the central region where Y* > 0.175 should fail 
structurally (Figure 2(b)). 

For the case of a friction angle of 35°, since uj* = 0.95 
and Y* is positive everywhere, the blue curve is always 
higher than, but close to, oj* (Figure 1(b)). Contrast to the 
25° friction angle case, at a spin rate of 0.9, since Y* is nega¬ 
tive at the central region, if Y is less 0.02, only the equatorial 
surface should fail (Figure 2(c)). At a spin rate of 1.15, on 
the other hand, and Y* is the highest at the central region 
(Figure 2(d)). For the case of a friction angle of 45°, w* for 
this case is 1.01. The results show the same features as the 


35° friction angle case; however, the region where F* > 0 
only appears at the equatorial edges. tJ*s for these cases are 
always lower than ui ~ 1.15. 

These cases describe that at high spin rates, a spher¬ 
ical body requires higher cohesive strength to avoid fail¬ 
ure of the central region. This corresponds to earlier stud¬ 
ies (Hirabayashi & Scheeres 2015; Hirabayashi et al. 2015). 
Plastic deformation should be compressive along the spin 
axis and tensile on the equatorial plane. However, if the spin 
rate is lower than , the equatorial surface is the most sen¬ 
sitive, which is one of the findings in this study. A variation 
of the friction angle causes a change of Y*- the higher the 
friction angle is, the smaller Y* is. As a result, oj* increases 
because a higher friction angle can prevent the body from 
failing at faster spin rates. For any cases, however, Y* of 
the equatorial surface {6 — 0°) becomes positive at slower 
spin rates than that of the central region; the failure regions 
move from the equatorial surface to the central region as the 
spin rate increases (Figure 1). This result indicates that the 
failure mode is not a strong function of a friction angle. It 
should be noted that this results from linear elasticity, which 
does not take into account plastic flow yet. Since plastic flow 
can make wider regions fail, the actual failure region should 
be wider than the current prediction (Chen & Han 1988). 
More precise propagation of plastic deformation in the body 
has to be discussed based on the flow rules, and we leave this 
discussion as our future work. 


4 NUMERICAL INVESTIGATION BY 

PLASTIC FEM 

We use a plastic FEM developed by Hirabayashi & Scheeres 
(2015) to evaluate the analytical results above. The code is 
developed by using commercial finite element program AN- 
SYS 15.03. In this model, elastic behavior is linear, while 
plastic behavior is perfectly plastic and follows an associ¬ 
ated flow rule, which defines that plastic deformation al¬ 
ways occurs in the normal direction to the yield surface. It 
is also assumed that plastic deformation is small. We cre¬ 
ate a 10-node tetrahedral FEM mesh for a spherical body. 
Six node displacements are fixed so that the rotational and 
translational motions are cancelled out. Similar to the ana¬ 
lytical model, we consider materials in the body to be ho¬ 
mogenous. The yield condition is the Drucker-Prager yield 
criterion, which is described in Equations (1) through (4). 

Figure 3 shows the results by the plastic FEM for the 
same cases as Figure 2. The contours indicate the stress 
ratio, a ratio of the current stress to the yield stress (Kohnke 
2009); if this ratio is 1, the region should experience plastic 
deformation, which is given in red. To induce the plastic 
deformation seen in this figure we choose the actual cohesive 
strength for each case. 

For the case of <j) = 25°, we choose Y = 0.02 at tu = 
0.9 to induce failure in both the central and surface regions 
(Figure 3(a)) and Y = 0.175 at w = 1.15 to produce that in 
the central region (Figure 3(b)). For the case of <j) = 35°, we 
input Y — 0.02 at to = 0.9 to trigger failure on the equatorial 
surface (Figure 2(c)), and Y = 0.15 at a; = 1.15 to make 
the central region fail (Figure 2(d)). Last, for the case of 
(j) = 45°, Y = 0.02 at Lj = 0.9, and Y = 0.125 at tj = 1.15. 
The case for tj = 0.9 does not have any failed regions in 


MNRAS 000, 1-10 (2002) 


Failure of spherical bodies due to rotation 5 




(a) (b) 



(c) 

Figure 1. Contour curves of Y*. The solid and dashed curves indicate Y* and Y*, respectively. Figures 1(a), 1(b) and 1(c) describe 
the cases oi (p = 25°, 35° and 45°, respectively. The red curves indicate the conditions of zero-cohesive strength for both cases. The blue 
lines give u} at which Y* = Y*. The color version is available online. 


the body, while that for oj = 1.15 causes the similar failure 
mode appearing in the central region. We confirm that the 
sizes of the failure regions by the plastic FEM are always 
wider than those by the analytical study. Regardless of this, 
all the results by both of the techniques are fairly consistent. 
We have to note that Figure 3 shows the deformed shapes 
after loading. 


5 DISRUPTION SCENARIO OF A 
SPHERICAL BODY 

The present analysis provided insights into evolution of a co¬ 
hesive, spherical body due to the YORP spin-up. Using the 
results obtained above, this section discusses possible evolu¬ 
tion and disruption scenarios of a spherical body by taking 
into account a two-layered-cohesive model by Hirabayashi 
et al. (2015). On the assumption that the bulk density is 
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(a) 



(c) 



(e) 



(b) 



(d) 



(f) 


Figure 2. Distribution of Y* across the quadrantal cross section along the x and z axes. The elastic solutions are used to provide this 
distribution. The left and right columns describe spin rates of 0.9 and 1.15, respectively. The first, second and third rows indicate friction 
angles of 25°, 35° and 45°, respectively. The color version is available online. 


homogeneous, their model consists of a spherical core which 
has higher cohesive strength and a surface shell which pos¬ 
sesses low cohesive strength. They concluded that if there is 
a strong internal core, a possible disruption mode is surface 
shedding. 

Figure 4 indicates evolution paths of a spherical body. 
Along every path (from the left to the right), the actual 
cohesion is considered to be constant, cu* and ut ~ 1.15 are 


used to characterize the stages of the evolution. A spherical 
body with homogeneous structure would have two pathways 
to its terminal state. If the actual cohesive strength is high 
enough for the body to support its structure at a slow spin 
rate, there are no failed regions in the body (Al). As the 
spin rate becomes higher than uic, the central region then 
reaches the yield first (A2). At this spin rate, although the 
failed region becomes wider, since the body can still sustain 
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(a) 



(c) 



(e) 


0.5 0.6 0.7 0.8 


ANSYS 



(f) 


0.9 0.95 0.99 1.005 


Figure 3. Failure regions on the cross section. The horizontal and vertical axes are identical to the x and z axes, respectively. The plastic 
FEM simulations are conduced to obtain these failure regions. The simulation conditions displayed here are the same as Figure 2. The 
color version is available online. 


itself, disruption does not occur yet. If the actual cohesive 
strength is not high enough, stage A1 can result in stage 
B2, which is discussed later. The body at stage A2 reaches 
its terminal state at a spin rate higher than to ~ 1.15 (A3). 
The actual cohesive strength is so high that the body can 
spin fast; however, once the structure fails structurally at 
a high spin rate, large components immediately fly away. 
This disruption mode is a breakup of the body into multiple 
components (Hirabayashi et al. 2015). 

If the actual cohesive strength is low at a low spin rate, 
the equatorial surface fails first (Bl). Stage B2 is a phase 
in which both surface and central regions fail at a spin rate 
higher than lj*, leading to propagation of plastic deforma¬ 


tion over the entire equatorial plane. This stage can come 
from two stages, Al and Bl. At this spin rate, since the equa¬ 
torial plane has already failed, the body becomes oblate due 
to vertical compression (Hirabayashi & Scheeres 2015). Once 
the spin rate is higher than uj ~ 1.15, the shape is more and 
more oblate (B3). The oblateness may depend on a friction 
angle because higher friction can prevent local elements from 
deforming inelastically (personal communication with Paul 
Sanchez, 2015). However, since plastic deformation highly 
depends on loading paths (Chen & Han 1988), if the body 
stays at critical conditions for a long time, the oblateness 
may grow critically even for high friction angle cases. Since 
the surface has failed already, a large amount of mass may 
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be ejected from the edges (Sanchez & Scheeres 2012). The 
prediction of this mode can be seen from Figures 2 and 3; if 
the actual cohesive strength is small, failure widely spreads 
out across the surface and central regions, leading to ver¬ 
tical compression (Hirabayashi & Scheeres 2015) and mass 
ejection. 

If the body has a strong core and weak surface, it has 
a different pathway to its terminal state (Hirabayashi et al. 
2015). At a spin rate of less than lj*, a failed region appears 
on the equatorial surface hrst (Cl). This mode should be 
identical to stage Bl. However, as the spin period becomes 
higher than u)c, since the central region is strong enough, 
the surface region suffers failure (C2). A considerable fail¬ 
ure mode at this stage is a landslide because the central 
region has high cohesive strength to resist failure, and the 
surface region should fail due to its low cohesive strength. 
The shape at stage C2 may look similar to that at stage B2. 
However, since the formation process is different, it may be 
possible to distinguish between stages B2 and C2 by observ¬ 
ing different surface morphology. The shape at stage B2 may 
not have the morphology created by landslides, while that 
at stage C2 may have. Once the spin rate is above oj ~ 1.15, 
only the surface mass is shed without large deformation of 
the original body (C3). This disruption can be distinguished 
from stage B3 because the amount of ejected mass and the 
hnal shape can be controlled by the size of the internal core. 
If the core is large, inelastic deformation is limited and only 
a limited amount of mass can be ejected from the surface. 


6 DISCUSSION 

The present analysis shows that the disruption modes of as¬ 
teroids may be dependent on the internal structure. On the 
assumption that reported active asteroids are spherical, the 
evolution diagram in Figure 4 shows that an asteroid hav¬ 
ing a breakup, such as P/2013 R3, should be structurally 
strong enough to fail at a high spin rate (A3). 62412 and 
P/2012 F5 (GIBBS) may be between stages B3 and C3 be¬ 
cause large components were not seen in their dust tails 
(Sheppard & Trujillo 2015; Drahus et al. 2015). Again, we 
emphasize that this is given only under the assumption that 
they are spherical at this point. In fact, they would be elon¬ 
gated rather than spherical. Further investigation is neces¬ 
sary to explore elongated cases. Another interesting aster¬ 
oid is 311P/PANSTARRS, formerly P/2013 P5, which has 
ejected a small amount of mass from the body (Jewitt et al. 
2013, 2015). This mode is consistent with stage C3. Also, 
the size of dust particles was only a few hundred microns. 
This gives rise to a possible scenario that the spin period of 
this body is not too fast, so ejected boulders may not have 
enough kinetic energy to escape from the main body. Note 
that our study does not explain the episodic event of this 
object. Again, further investigation is necessary to explain 
this unique event. 

It is interesting to consider possible formation scenar¬ 
ios for a quasi-spherical object with equatorial ridges. Either 
stage B2 or stage C2 may represent the observed shapes such 
as 1999 KW4 Alpha (Ostro et al. 2006), 1950 DA (Busch 
et al. 2007) and 2008 EV5 (Busch et al. 2011). To have 
these stages, they should go through one of the following 
processes: Al —> B2, Bl —> B2 and Cl —^ C2. To go through 


process Bl —>■ B2, a homogeneous body has to have small 
cohesive strength. Based on Figure 2(c), for a friction an¬ 
gle of 35°, a spherical body should have an actual cohesive 
strength of less than 0.02 at a spin rate of 0.9. For these 
asteroids, we find that the actual cohesive strengths of 1999 
KW4 Alpha, 1950 DA and 2008 EV5 should be less than 
8 Pa, 5 Pa and 0.5 Pa, respectively. For 2008 EV5, we as¬ 
sume its bulk density to be 2000 kg/m®. Note that these 
values are quite small compared to the reported cohesive 
strength of 1950 DA, which is > 75 — 85 Pa (Hirabayashi & 
Scheeres 2015), and that of P/2013 R3, which is 40 — 210 
Pa (Hirabayashi et al. 2014). 

We also investigate w/ for these asteroids. The central 
region begins to deform at this spin state and becomes more 
oblate as the spin rate increases. Assuming a friction angle 
to be 35°, we obtain this boundary as ~ 3 hours for 1999 
KW4 Alpha, 1950 DA and 2008 EV5. Since the spin periods 
of 1999 KW4 Alpha and 1950 DA are 2.7 hours and 2.1 
hours (Ostro et al. 2006; Busch et al. 2007), respectively, 
these asteroids can become more oblate due to failure of 
their central regions at the current spin states. On the other 
hand, the spin state of 2008 EV5, which is 3.8 hours (Busch 
et al. 2011), is below this lower threshold, implying that the 
body may structurally relax. However, because of its shape, 
it might have been spinning at a spin period of shorter than 
3 hours. 

Finally, we mention the earlier studies about the shape 
evolution of a cohesionless, spherical body. Harris et al. 
(2009) developed a numerical model that computes an ideal 
shape of a spherical body by considering the angle of repose. 
Extending their study analytically, Scheeres (2015) found 
that in the two-dimensional model, once failure occurs on 
the surface, the material that comes off of the mid-latitudes 
should leave behind material at the angle of repose for the 
regolith. This is different from our results that showed that 
the equatorial surface should fail hrst at a slow spin rate. 
The difference between our work and his may come from 
our computational settings. Specihcally, our boundary con¬ 
dition may be different from his. Our surface condition is 
based on a zero-traction condition (Hirabayashi et al. 2015), 
while he dehned the surface stress based on the force bal¬ 
ance on the surface. These are not necessarily the same. We 
will investigate the connection between our work and his in 
the future. 


7 CONCLUSION 

This paper explored how failed regions in a spherical body 
evolve as the spin rate changes and deduced its termi¬ 
nal state based on our results and the earlier works by 
Hirabayashi & Scheeres (2015) and by Hirabayashi et al. 
(2015). The main contributions are as follows: 

• Prediction of the failed regions by the elastic model is 
consistent with the results by a plastic FEM. 

• The surface region at the equator needs higher cohesive 
strength to avoid failure at a low spin rate, while the central 
region must have much higher cohesive strength as the spin 
rate increases. 

• There is a spin rate condition, given as w/, at which the 
central region needs cohesive strength to avoid failure. If the 
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Figure 4. Evolution and disruption scenarios of a spherical body. The blue and yellow markers show elastic components, and the orange 
markers describe failed ones. The blue components are structurally weaker than the yellow components. As the spin rate increases, the 
disruption mode evolves differently due to the internal structure. The color version is available online. 


spin rate is higher than this condition, a spherical body may 
become oblate. 

• OJ* increases as the friction angle becomes higher; the 
friction angle changes the failure conditions. However, even 
at different friction angles, the failed regions consistently 
move from the equatorial surface to the central region as the 
spin rate increases; the failure modes are not strong func¬ 
tions of it. 

• Three possible disruptions were discussed in this study. 
For a homogeneous body, a strong structure allows the body 
to remain at a higher spin rate and then causes it to have a 
breakup into multiple components. A weak structure leads 
to failure on the equatorial surface at a slow spin rate, and 
then has failure propagation to the central region. As a re¬ 
sult, at its terminal state, the body becomes oblate. For a 
body with a strong core, a possible failure mode is only sur¬ 
face shedding. 

• The formation of the equatorial ridges on a spherical 
body depends on cohesive strength. If the surface region 
is very weak, the ridge should be made only by downslope 
movement, which mainly occurs at lower latitudes; other¬ 
wise, the main driver for it is inelastic deformation along 
the spin axis. 
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